No many-scallop theorem: Collective locomotion of reciprocal swimmers 
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O 

' To achieve propulsion at low Reynolds number, a swimmer must deform in a way that is not 

invariant under time-reversal symmetry; this result is known as the scallop theorem. We show 
here that there is no many-scallop theorem. We demonstrate that two active particles undergoing 
reciprocal deformations can swim collectively; moreover, polar particles also experience effective 
1/^ , long-range interactions. These results are derived for a minimal dimers model, and generalized to 

more complex geometries on the basis of symmetry and scaling arguments. We explain how such 
' cooperative locomotion can be realized experimentally by shaking a collection of soft particles with 

^ ' a homogeneous external field. 

■ PACS numbers: 87.19.ru, 47.15. G-, 47.63.-b, 62.25.-g 



' O ' Microorganisms rely on their ability to swim in order to achieve a variety of biological tasks, including sensing, 
Q [ targeting or feeding [l| . Well-studied examples include bacteria, spermatozoa and ciliated protozoa . Given the 

^ ' recent advances in the construction of complex colloidal assembly P, Q and in the coupling of biological machines to 
artificial microstructures [6|] , man-made functional micro-swimmers are expected to catch up with real microorganisms 
^ , ItI. However, a fundamental challenge in designing artificial micro-swimmers lies in the constraints of the so-called 

cn , " n n 

0^ ■ scallop theorem [8|, |9|. Since at small scales, or at low Reynolds number, the Stokes flow equations are linear and 
time reversible, swimming can only be achieved by a sequence of shape deformations invariant under time-reversal, 
or reciprocal; the prototypical reciprocal movement is that of a scallop which opens and closes its shell (Fig.[Tl'V, left). 
From a design perspective, the scallop theorem implies that, in order to move on average, a low-Reynolds number 
swimmer must have at least two internal degrees of freedom undergoing a non-reciprocal cycle (Fig. [TJ'V, right). 

_ ^ ' In this paper, we show that although an isolated reciprocal active particle cannot move at low Reynolds number, 
' two or more are able to swim collectively; in other words, there is no many-scallop theorem. We demonstrate that 
active particles modulate the relative distance by exploiting hydrodynamic coupling, thereby inducing on each other 
velocity fields which do not average out to zero and allowing for collective motion to occur. We also show that a set 
of non-identical reciprocal active particles experience long-range effective interactions which can either be attractive 
or repulsive depending on the particles geometry. These results are demonstrated rigorously for a pair of force- 
free reciprocal dimers, and extended to more general geometries on the basis of symmetry principles and scaling 
arguments. Experimentally, we show that simple elastic particles shaken by a homogeneous oscillating external field 



* Email: elaugaSucsd.edu 

t Email: denis .bartoloOespci . f r 



Reciprocal Non-reciprocal 

(A) (y) =0 (^) 7^ 




fa fa fb fb 



(B) 

Gil 



d{t) 



Lit) hit) 



FIG. 1: (A) A body deforming its shape in a reciprocal fashion, such as a scallop, cannot move on average at low Reynolds 
numbers (left, (V) — 0), whereas non-reciprocal deformation leads to net propulsion (right, (V) 7^ 0) ; (B) Two force- free dimers 
interacting hydrodynamically. The dimers are composed of two solid spheres, of radii {ai, a2} and {fei, 62}, have lengths £a{t) 
and £bit), and are separated by the distance d{t); the forces on each sphere in a dimer are equal and opposite. 

can be exploited to obtain collective locomotion. 

We consider a collection of prototypical active particles, force-free dimers, for which we neglect the flow disturbance 
created by the links joining them (see Fig. [T}3). The i'^ dimer is composed of two spheres, ii and i2, separated by 
the time-varying distance £i{t) = Xi^ — x^j, where Xa denotes the position of sphere a = ji,«2 along the a: axis. The 
linearity of Stokes equation implies that the motion of each sphere in a dimer is linearly related to the forces, fp, 



acting on the spheres from all dimers 10 1 



Xa=^Haj3ff!, (1) 

where the {Hap} are the hydrodynamic mobilities of the spheres {(3 = ji, j2)- We assume that each dimer is force-free, 
fii + /i2 ~ 0' ^^^d define fi = fi^ = —fi^ (Fig. [TJB). The center of the i**^ dimer is denoted Xi = {xi^ + Xi^)/2. 

Single dimer. We first consider the case of an isolated active dimer, labeled a. In that case, we have Xa = H{ia)fa 
and the dimer elongation satisfies £a — fJ-a{ia)fa- Hence, we have Xa = H{£a)£a/ fJ-ai^a), which is an exact derivative 
and therefore averages to zero over time, (xa) = 0, for any periodic sequence of forces or deformations imposed to the 
dimer. This is the scallop theorem. 

Two dimers. We now consider the case of two dimers (a and b). We demonstrate that, although the dimers are not 
able to swim when alone, hydrodynamic interactions through the viscous fluid enable collective motion. We restrain 
our analysis to widely separated dimers in an unbounded fluid; if d is their relative distance, we consider the limit 
£i ^ d in which case the mobilities in Eq. ([1]) are given by the Oseen tensor [iq|: Ha[3 ~ l/Ar/nd = H{d), if a ^ /3, 
and Haa = l/^rjiraa = H{a) otherwise (77 is the fluid viscosity, and the radius of the sphere). Performing a Taylor 
expansion for the dynamics of the four spheres, Eq. [l] and keeping only terms of order 1/d^ lead to a system of 
four nonlinear differential equations, two for the positional degrees of freedom, Xi, and two for the internal degrees of 
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freedom, ii, as given by 



Xa = Mafa - [dxH{xb - Xa)]hfb, 

Xb = Mbfb + [dxH{xb - Xa)]tafa, 

£a = ^J■aiia)fa + ^J.ab(aibfb, 

4 = tJ'b{ib)fb + tJ'abidbfa, 



(2a) 
(2b) 
(2c) 
(2d) 



where we have defined Mi = [H[ii) — H{i2)\ /2, ^ab = dxxH{d) and ^^(4) — 2H{ii) — H[ii) — H{i2), mobUity 
coefficients associated with the position and the elongation of each dimer respectively. We further restrain our 
analysis to small amplitude reciprocal motion of each dimer around the time average length, ii = {£i) (in all that 
follows brackets stand for time average). More precisely, if we write £i — £i + 0(e, e^) and fi ~ 0(e,e^), we keep 
terms up to O(e^) in Eqs. (Pa|l - (l2dp . We can then compute the average collective and relative swimming speeds of 
the dimers, (V) = {xb + Xa)/2, and AV = {xb — Xa) respectively. To do so, we distinguish the two ways in which the 
dimers can be physically actuated. 

Force-driven motion. We first consider the case where the internal forces, /i, are specified. This is analogous 
to biological swimmers possessing force-generating units (the axoneme for eukaryotic cells the rotary motor for 
bacteria such as E. coli ,12]). In that case, we assume the internal force to be known, 0(e), and time-periodic. The 
force-displacement relation, Eq. ((2c|) - ([2d|) . can then be linearized, and after some straightforward algebra, we obtain 
for a pair of dimers in an unbounded fluid 
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Generalization for an arbitrary response function H, including for example boundary effects, is given in 121]. The 
results of Eq. ([3]) show that, generically, the scallop theorem breaks down for two active particles interacting hydro- 
dynamically: taken individually, these particles cannot move, but when interacting through the fluid, they display 
collective motion {{V) ^ 0), and experience an effective long-range interaction {AV ^ 0), both of which decays in 
space as The direction and sign of the collective and relative speeds depend on the dimers' geometry and 

actuation; as the simplest example, if we consider sinusoidal forcing of the form fi{t) ~ fiCOs{ujt -\- we have 
if a J fb) = /a/6sin((/)b — (j)a)/2Lj, and locomotion occurs if the two particles are actuated with phase differences, a 
direct consequence of the time reversal symmetry of the Stokes equation. 

Physically, locomotion of the pair of dimers occurs because their relative distance is oscillating in time, and therefore 
the flow fields seen by each dimer does not average to zero. In fact, if the dimers are rigidly connected (i.e. if d(t) is 
kept constant), we would obtain {V) = at order l/d'^, and therefore fluid- mediated forces are the crucial ingredient 
leading to collective locomotion. 

Displacement- driven motion. We now assume the sequence of deformation of each dimer, 4(i), to be specified and 
time periodic, 4=4 + S£i{t), with 5£i{t) — 0{e). This is the relevant limit for (robotic) man-made micro-swimmers. 
We now have to invert the force-displacement relationship, Eq. (|^ - (|2dp . and after some tedious but straightforward 




FIG. 2: Dynamics of identical dimers interacting hydrodynamically. Units are chosen so that Slait) — cos{t), S£i,{t) — sm{t). A: 
mean displacement, (J V dt), for two mirror- images polar dimers (ai 7^ 02, solid line); mean displacement for two mirror-images 
apolar dimers (ai — a2, dash-dotted line); mean distance, J AV dt, between the two dimers (polar and apolar cases, dashed 
line). B: mean displacement (solid line) and mean separation distance (dashed line) for two identical polar dimers pointing 
in the same direction. C: same as [2)3 but for dimers pointing in the opposite direction; note the change in the sign of the 
dimer-dimer effective interaction. D: comparison between the swimming speed of two active dimers and a single swimmer 
made of two rigidly connected dimers (constant d). In all figures, the lengths £a and it oscillate with the same amplitude and 
frequencies, and a relative phase of 7r/2. In the chosen units, £a = it ~ 10, d\t=o ~ 100, ai — 3, a2 = 3/2 for the polar dimers 
and tti = 02 = 2 in the apolar case. 
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algebra we obtain locomotion at speed 

{V) = 
AV = 

where we have defined i = 1112/ (ii +^2) (i — a, b). The general expressions for an arbitrary kernel H are given in 22 1. 
Similarly to force-driven motion, the scallop theorem does not hold for a pair of reciprocal bodies, and locomotion 
arises if there is a non-zero phase difference between the deformation of each particle. If the dimers were rigidly 
connected, we would obtain in this case a Ol\/(P) swimming velocity, but with a sign opposite to that arising from 
hydrodynamic interactions, Eq. (Pa|) [23! . 

Identical swimmers. In the particular case where the two dimers are identical (a,; — bi, i — 1, 2, and £a — = ^)) 
the result of Eqs. ([5a|) - (|^a|l cancel out. One needs to go to higher order in the asymptotic expansions to obtain the 
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swimming kinematics, which is readily done, and locomotion is obtained at order l/d^, with velocities 

- lQnW~ad^ ' = ' 

for force- and displacement-driven motion respectively. 

In Fig. [2]A_-C we highlight the main features of the collective dynamics of a pair of identical dimers in the case of 
displacement-driven motion, by integrating numerically Eqs. ()4bp and ([5]); results for force-driven motion are similar. 
We first show in Fig. [2]A that, on average, mirror-symmetric swimmers remain at the same relative distance and 
swim collectively with a constant velocity. When the swimmers do not display mirror-image symmetry, we show in 
Fig. [2}3-C that they undergo a repulsive or an attractive effective interaction depending on their relative orientation. 
The dimers are separating/approaching at rate d ~ il/d"^, and therefore d ^ {to ±t)^^^; accordingly, the mean 
speed decreases as (V) ~ (to ± t)~^, and the distance traveled by the swimmers can be arbitrarily (logarithmically) 
large, J {V) dt ^ log(to ± t) 24]. Note that in the case of relative attraction, a proper description of the near-contact 
hydrodynamics would regularize the finite-time singularity displayed in Fig.[2j3; in an experiment, short-range surface 
forces (such as van der Waals) would lead to self-assembly of the two dimers into a four-spheres non-reciprocal 
swimmer. 

Generalization. The results above demonstrate the emergence of collective motion and of long-range interactions 
between two dimers embedded in a viscous fluid, and can be generalized in a number of ways. 

First, our calculations were made under the assumption of an infinite fluid, but the results remain valid in other 
geometrical settings, such as a Hele-shaw cell, or in confined geometries typical of microfluidic applications, by changing 
appropriately the Oseen-like kernel H{d). In addition, our dimers are made of two connected spheres, but other solid 
bodies with l/a ^ 1 can be considered by simply modifying the values of H{a) [21,22]. 

Furthermore, Eqs. ([3]) and ([4]) also apply for stochastic fluctuations of the length of each dimer, and any correlation 
between the noisy shaking sources acting on each dimer is seen to lead to non-zero swimming velocities and effective 
interactions between the two active objects. 

Theoretically, the 1/d'^ spatial decay of the velocities with the intcr-dimer distance in Eqs. ([3|) and ([4]) could have 
been anticipated since the locomotion arises from the rectification of interacting force dipoles decaying as at 
long distance ^13|; this scaling argument also hold for any pair of active particles having a single deformation degree 
of freedom aligned with the a;-axis. In the small deformations limit, their collective and relative swimming velocities 
is expected to be of the generic form: 



{V) ~ (/,//,)^^, AV^{Mh)J2^, (6) 

n>3 n>3 

{V) ^ {SiM AV^ {6i,SL) E (7) 

n>3 n>3 

for force- and displacement-driven motion respectively, and where /3n, 7n, f^n} depend solely on the shape of the 
particles [25^. Note from Eqs. (|3bp - (|4b[) that for two mirror-image particles, 73 = (^3 = 0; this result is actually true 
at all orders, and two swimmer with mirror-image symmetry verify AV = 0. Indeed the flow and pressure fields 
induced by the beating of two mirror-images dimer are invariant under the combination of time-reversal and mirror 
symmetries. Consequently, AV, which has the symmetry of a velocity gradient, transforms into —Ay, and therefore 
AV = 0. Similarly, we have ~ (3^ ^ Q for identical dimers; this is a consequence of the two (rectified) dipoles 
having opposite contributions on each dimer, and is generally true for a pair of identical particles. 
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Beyond the simple one-dimensional models considered here, all types of reciprocal motion (including those with 
three-dimensional shape deformation) are also expected to display collective motion induced by hydrodynamic inter- 
actions. And although we have emphasized locomotion in this paper, our results can be extended to the dual problem 
of pumping fluid by anchored bodies. 

Soft swimmers. We now turn to the experimental realization of these ideas. Actuating a collection of active 
particles with out-of-phase conformational changes is difficult, and it would be preferable to devise a framework where 
an homogeneous external field could produce directed locomotion. We show below that this is possible if the particle 
are soft: indeed, since soft particles in a viscous fluid posses relaxation times, different particles will naturally react to 
an external shaking with phase differences, and locomotion will ensue. Practical examples of soft-particle actuation 
include bending of magnetic fllaments by magnetic fields , tem pera ture or light actuation of liquid-crystal elastomers 
Q, [3 , or self-sustained chemical reactions for swelling of gels jl6| . 

We thus consider a pair of force-free apolar and identical dimers (radius a, average length £) subject to a homo- 
geneous external shaking. We write fi = /^^ako _^ f^'^^^^, for i = a,b, where the force f^^'^^" is externally produced 
and the same for each dimer (homogeneous forcing). The force f'^'^^^^ is the internal (elastic) response of the dimer, 
and write — kiSii where ki is the dimer stiffness; in that case, the intrinsic dynamics of each active particle is 

characterized by a relaxation time scale — Snria/ki. Assuming a monochromatic shaking f^^^^'^ — /o cos(ajt) for 
simplicity, integration of Eqs. (I2cl) - (j2dp leads to collective motion with speed 



(V) 



(1+C.V2)(1+..V2) 



(8) 



'b . 

and Ay = by symmetry. We see that, under homogeneous forcing, the only condition necessary to obtain locomotion 
is Tq ^ Tfc. Locomotion always occurs in the direction of decreasing relaxation time, i.e. the stiff dimer is pulling the 
soft one; optimal locomotion occurs when the system is actuated with frequency lo ^ ija'Th)^^^'^ and when the ratio 
of relaxation times is large. 

A final relevant example is the case where one dimer is purely passive, i.e. its elongation is not coupled to the 
external field. Specifically, we assume dimer a to be actively shaken, with f^^^^° = fo cos(ct't), while dimer b is passive, 
and ff^^^'^ ~ 0. In that case, the beating of the passive dimer arises from the hydrodynamic interactions with the 
active dimer, and its amplitude is a function of the distance d. Collective motion arises with speed 



327r2?72d7' 



(9) 



.(l+c.V2)(l+a;V2). 

and the active particle is seen to "pull" the passive one. 

Perspective. We have shown in this paper that there is no many-scallop theorem: two bodies with reciprocal de- 
formation can exploit hydrodynamic interactions to obtain collective locomotion and effective long-range interactions. 
On the basis of our two-body investigation, the generalization to a large number of reciprocal soft active particles 
interacting through a fluid is expected to display both spatial and spatio-temporal spontaneous symmetry-breaking 
(coarsening, ordering and directed large-scale motions). Similar cooperative effects for coupled active (but non-self 
propelled) particles have already been demonstrated in the very different context of molecular motors 17[, active 
membranes dynamics [3| and vibrated nematic grains [foj l. 

This work was supported in part by the NSF (grant CTS-0624830 to EL), by ESPCI and by Denis Diderot University 
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